Cooper pairs as resonances 
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Abstract 

Using the Bethe-Salpeter (BS) equation, Cooper pair- 
ing can be generalized to include contributions from holes 
as well as particles from the ground state of either an ideal 
Fermi gas (IFG) or of a BCS many-fermion state. The 
BCS model interfermion interaction is employed through- 
out. In contrast to the better-known original Cooper pair 
problem for either two particles or two holes, the general- 
ized Cooper equation in the IFG case has no real-energy 
solutions. Rather, it possesses two complex-conjugate so- 
lutions with purely imaginary energies. This implies that 
the IFG ground state is unstable when an attractive inter- 
action is switched on. However, solving the BS equation 
for the BCS ground state reveals two types of real solu- 
tions: one describing moving (i.e., having nonzero total, 
or center-of-mass, momenta) Cooper pairs as resonances 
(or bound composite particles with a finite lifetime), and 
another exhibiting superconducting collective excitations 
analogous to Anderson-Bogoliubov-Higgs RPA modes. A 
Bose-Einstein-condensation-based picture of superconduc- 
tivity is addressed. 
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1. Introduction 

The original Cooper pair (CP) equation is a 
two-electron Schrodinger equation in momentum rep- 
resentation with a given two-body interaction (having 
some attraction) but includes ad hoc restrictions on 
the magnitudes of both electron wavevectors ki, k2, 
namely ki > k^, k2 > k^, where kj, is the electron 
Fermi wavenumber for an ideal nonrelativistic many- 
electron system. One then seeks the energy eigenval- 
ues of a CP bound state and its corresponding wave- 



function. 

Since there is no rigorous derivation of the original 
CP equation, several authors (^-||] reformulated the 
complete CP problem without neglecting holes, using 
the mathematically exact BS equation approach H 
applied to the system, in search of two-particle bound 
states in the presence of other system electrons. Such 
a treatment allows generalizing several approaches 
in superconductivity theory, including the BCS, the 
BCS-Bose crossover, and the Bose-Einstein condensa- 
tion pictures. But if the BS equation is based merely 
on the IFG ground state one obtains purely imaginary 
solutions, suggesting that the ground state is unsta- 
ble in the presence of attractive interactions of some 
kind. This is an instability of the IFG ground state 
with respect to the creation of two-particle (2p-) or 
two-hole (2h-) resonant states, a situation analogous 
to the classical problem of hydrodynamic instability 

i- 

2. BS equation based on BCS ground state 

Consider, however, the generalized two-component, 
two-electron BS equation based not on the IFG ground 
state [||-|] but rather on the BCS ground state. We 
introduce the Bogoliubov-Valatin u, v transformation 
of electron Fermi operators ak,cr to new Fermi opera- 
tors Qfk, cr, namely 

aic.o- = Wfc "k.o- + 2(T Vfc a+jj. (1) 

where cr = ±i is the spin projection for an electron 
state. The coefficients Uk,Vk are real and depend on 
k. The many-electron system hamiltonian is H = 
Hq + Hint where 
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where Ck = h^k^ /2m is the kinetic energy of an elec- 
tron of mass m; Ep = h'^kp/2m is the Fermi en- 
ergy; L is the system size; v(q) the Fourier trans- 
form of the two-electron interaction potential; and 
the last sum is restricted by momentum conservation 
k[ + k'2 = ki + ko. This leads to the well-known BCS 
hamiltonian 

Hbcs ^Uo + J2^ ^E{k)a+ ^ak, (3) 

with Uq a generalized BCS ground-state energy 
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E{k) = y/A^{k)+B^{k). (4) 



Here B{k) plays the role of the original BCS energy 
gap A(fc), and 

A{k) = ek-EF + 

+ 2L-M0)Y.<' - L-'Y.''{\k - k'|)4 , 

k' k' 

i3(fc)EEL-3EK|k-k'|K.i;,,. (5) 
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For the new interaction hamiltonian we obtain 
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x2cri Q:-ki,-o-iQ:ki(Ti2CT2 a-k^,-(T2Q!k20-2 ■ 
where A:') = Uk Uk' — Vk Vk' and M{k, k') = 

Uk Vk' + Uk' Vk . 

To obtain the BS equation based on the BCS ground 
state consider the Feynman diagrams of perturbation 
theory based on this ground state, where (y|) is the new 
unperturbed hamiltonian Hq. We now have the usual 
arrowed electron lines labeled by k, E, a to which we 
associate the BCS unperturbed Green's function 



go(k,E,a)^- ^ 



i ~E + E{k) - ie' 



(6) 



where E{k) is given by (jj). There exist four-line- 
end double vertices of six different kinds (see Fig. 
1) where the interfermion interaction is denoted by 
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dashed lines. To a double vertex type (a) of Fig. 1, 
with two outgoing line ends with indices (k'j^, E'l, ai) 
and (kj, i?2, (^2) along with two incoming line ends 
with indices (ki, Ei, ai) and (k2, £'2, 0-2), we attach 
the factor 

i.(|ki - k'll) L(ki, k'l) i(k2, k'2). 

To a double vertex of type (b), with two outgoing line 
ends (k'^j Ei' , (Ti)and (kj, E2' , CF2) along with two 
outgoing line ends (— ki, — — cti) and (— k2, — -E2, 
— 1T2), as well as to a double vertex of type (c) with 
two incoming line ends with indices (— k'j^, —Ei , — (Ti) 
and (— k2, —E2' , —(72) along with two incoming line 
ends with indices ki, Ei, ai) and (k2, E2, 172), we 
attach the factor 

-L-^i^dki -k'll) 4(71(72 M(ki, k;)M(k2, k'2). 
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FIG. 1. Six different types of vertices with correspond- 
ing topological automorphism factors g. 

3. Coupled BS equations 

Because of different kinds of vertices in Fig. 1 we 
now have a system of two coupled BS equations. Fig- 
ure 2 shows their diagrammatic representation in the 
ladder approximation of the two-electron BS equation 
for the BCS ground state. Depicted are both the two- 
electron V'+lk-E; K£/<-) and two-hole K£/j-) 
bound-state functions. Here K = ki -I- k2 is the total 
(or center-of-mass) wave-vector and £k is the total 
energy of the two electrons referred to 2Ep. 

Using the diagramatic rules just described, the two- 
component BS equations for the bound state is 

^+(hE)=-— / dE K|k-k'|)x 



X L{K/2 + k, K/2 + k') L(K/2 - k, K/2 - k')x 
X {i/hfGo (K/2 -I- k,£A72 + E)x 
xGo (K/2 - k,£K/2 - E) V+(k'£;')- 

-— / dE K|k-k'|)x 

X Af (K/2 + k, K/2 + k') M(K/2 - k, K/2 - k') x 

X {t/hfOo (K/2 + k,£K/2 + E)x 
xgoiK/2-k,£K/2-E)^P^{kE'); (7) 

^_lkE) = -— j dE L-^Y.^,v{\k-k!\)x 

X L{K/2 + k, K/2 + k') L(K/2 - k, K/2 - k')x 
X {t/hfgo (-K/2 - k, -£k/2 -E)x 
xgn{-'K/2 + -k,-EK/2 + E) i!^{kE')- 

-^J_^dE L-^Y..,.ilk-k'l)x 

xM(K/2 + k, K/2 + k') Af (K/2 - k,K/2 - k')x 

X (i/hfga i-K/2 - k, -£k/2 - E) x 

X go(-K/2 + k,£A72 + i^) ^+(k'£;'). (8) 

It can be shown that these equations coincide exactly 
with those used in the description of collective excita- 
tions in the BCS-Bogoliubov microscopic theory 
of superconductivity. 
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FIG. 2. Diagrammatic representation of the two-com- 
ponent coupled BS equations. 
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4. CP solutions of BS equations 

We now employ the BCS model interaction 

K|k-k'|)^-(4/A:'')F?7(fc)r7(fc'), 



(9) 



where V > 0, and ri{k) = 1 when kp — ki) < k < 
kp + kjj and = otherwise. Here kjj = mujjj/hkF 
with ujd the Debye frequency, if Hcod << Ep. 

The detailed solution for the system of the two- 
component coupled BS equations of Fig. 2 is too cum- 
bersome to present here, but it can be shown that they 
yield two types of independent solutions. The first so- 
lution of (0) and (H) with is just 



Vk 
"2 
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xi;(|K/2 - k|) [u(|K/2 - k|) w(|K/2 - 

- u(|K/2 + k|)i;(|K/2-k|)] > 

^(|K/2 + k|)+^(|K/2-k|) 
-fi + [i?(|K/2 + k|) + i?(|K/2-k|)P 
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' being the angle between k 
k| = yW+KktTlOjl and 



Here t = cos 6*, 
and K while |K/2 
|K/2 - k| = - Kkt + KlJA. 

As we see the numerator in (|l^) vanishes for X = 0, 
so the denominator must vanish as if ^ as well. 
This gives £k = ±2 A -h 0{K). To find the asymp- 
totic solution for both small coupling and K we may 
neglect terms with A in u, w in the integrand. Thus 
we put u{k) ~ e{kF -k) = Opik), v{k) ~ e{k - kp) = 
Oaik), where 9{k) is the usual Heaviside step function. 
Therefore, 

u(|K/2-|-k|)w(|K/2-k|) x 
X [u(|K/2-k|)i)(|K/2 + k|)- 
-7.(|K/2-fk|)^;(|K/2-k|)]c. 

-^f(|K/2 + k|)(?G(|K/2-k|). 
Thus we can rewrite (|l0| ) as 

f / dk dtx 

xeF{\K/2 + k\)9G{\K/2-k\)x 
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£k + E (|K/2 + li\) + E (|K/2 - k|; 



= 1. (11) 



In Fig. 3 we show as shaded the region of integration 
over k and t = cos 6. The integral over the polar angle 
9 is restricted to from to 7r/2 since the integrand in 
( |lO|) vanishes for t < 0. We first integrate over k from 
kmin to fcmax which are solutions of |K/2 + k| = kp, 
and |K/2 - k| = kp, or 



^Kt/2 + K^/A + kj,, 



kp . 



The restrictions kp — kjj < k < kp + kjj are not 
illustrated in Fig. 3 as they are unimportant for us. 
As if ^ we then have 



-Kt/2 + 0{K'^), 
Kt/2 + 0{K^), 
E (|K/2 + k\)+E (|K/2 - k|) ~ 2 A + 0{K^) 



^min 
k 



where we redefined the integration variable k ^ kp 
+xK^ with X the new variable. 
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FIG. 3. Integration region in 

Consequently in seeking a solution of Eq. (^ij), for 
example of the form 

£k ^2A + chK + 0{K^), 

where c is a constant, we obtain for small A and K 

' dt — ~ 1. 



An'^h Jo 



kF~Kt/2 



cK 



Thus c = V k'^/S-n'^h. Finally one gets 



£k ^2A+ -vphK + O(K^) 
4 



(12) 
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together with its symmetric solution 



superconductivity ^0 
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Sk = -vphK + 0{K^), (13) 

where A = N{0)V and A^(0) = mkp /l-rr'^h^ ■ These 
two solutions describe moving 2p-CPs and 2/i-CPs, re- 
spectively. A linear -'m-K behavior of the moving pair 
binding energ y w as obtained for the original Cooper- 
pair problem [[ll| but it was independent of the inter- 
action coupling strength as it excluded 2h-CP contri- 
butions. More significantly, the binding energy there 
was negative as it refers to an infinite- lifetime compos- 
ite particle, while in ( p^ ) it is positive as it describes 
a resonance in the continuum with a finite lifetime as 
evidenced by an imaginary contribution [T^ ] appear- 
ing in higher order terms in K . 

5. ABH-like mode solution of BS equations 

The second solution of the coupled BS equation (|^) 
and (^) follows from 

^ / dk / dtt.(|K/2 + k|)x 

27'' JkF-ko J-1 

xv{\K/2 - k|)[u(|K/2 + k|) u(|K/2 - k|)-H 



6. Conclusion 

We have presented a new many-fermion formalism 
based on a BS equation applied to the full BCS ground 
state which does not neglect the presence of holes. 
This leads in the ladder approximation to two types 
of solutions: the first referring to simple moving CPs 
consisting of two-electron (12) or two-hole (13) reso- 
nances. In addition, this formalism naturally provides 
a second solution (15) of an entirely different physi- 
cal nature which is analogous to the ABH excitation 
mode. Bose-Einstein condensation can occur with the 
first type of objects but not with the second. 
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Proceeding as before, in the limit K ^ one now 
finds 

£k = {VF hK/V3)[l + J|^e-2A + ...]+ o{K'), 

(15) 

which is similar to the Anderson-Bogoliubov-Higgs 
(ABH) RPA excitation mode in the BCS theory of 
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